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The Integrate-and-Dump Filter (IDF) is used as a matched filter for the detection of 
signals in additive white Gaussian noise . In this article , the performance of the digital 
integrate-and-dump filter is evaluated. The case considered is when symbol times are 
known and the sampling clock is free running at a constant rate, i.e., the sampling clock 
is not phase locked to the symbol clock . Degradations in the output signal-to-noise ratio 
of the digital implementation due to sampling rate, sampling offset, and finite bandwidth, 
resulting from the anti-aliasing low-pass prefilter, are computed and compared with those 
of the analog counterpart. It is shown that the digital IDF performs within 0.6 dB of the 
ideal analog IDF whenever the prefilter bandwidth exceeds four times the symbol rate 
and when sampling is performed at the Nyquist rate . The loss can be reduced to 0.3 dB 
by doubling the sampling rate, where 0.2 dB loss results from finite bandwidth and 
0.1 dB results from the digital IDF. 


I. Introduction 

An Integrate-and-Dump Filter (IDF) is the ideal matched 
filter for coherent detection of rectangular pulse shape signals 
corrupted by additive white Gaussian noise (AWGN). A digital 
implementation of the IDF has numerous advantages over its 
analog counterpart, such as the ability to dump instanta- 
neously with no overshoot, freedom of drift from the quies- 
cent operating point, and the use of advanced digital inte- 
grated circuits (ICs) to perform multiplication and accumula- 
tion with greater accuracy than the analog counterparts. 

A. Current Problem 

The digital implementation of the IDF requires the input 
waveform to be sampled. This requires that the input signal be 


filtered prior to sampling to eliminate aliasing. Throughout 
this article, the anti-aliasing low-pass prefilter is also referred 
to simply as the filter. Filtering of the observed signal results 
in the transmitted signal pulses becoming band limited and 
also causes Inter-Symbol Interference (ISI). The finite band- 
width of the observed signal and ISI both degrade the perfor- 
mance of the receiver. In order to evaluate the performance 
versus the filter bandwidth and sampling frequency, the analog 
IDF depicted in Fig. 1(a) is compared with the performance of 
the digital IDF shown in Fig. 1(b). 

In the limit, as the prefilter low-pass bandwidth W Hz and, 
as a consequence, the sampling rate approach infinity, the per- 
formance of the digital IDF converges to that of the analog 
IDF. However, this is not practical from an implementation 
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point of view, since large bandwidths require higher sampling 
rates to satisfy the Nyquist criterion. In this article, we find 
the performance degradation which results from filtering the 
received signal and from using only a small finite number of 
I samples. 

| Furthermore, we also superimpose the effect of “offset 
sampling” in the performance of the digital IDF. When sam- 
pling a signal of finite duration T sec with a sampling period of 
| T s sec, the first sample of the signal may occur anywhere in 
i the time interval 0 < t < T y This consideration for the digital 
IDF has not been studied previously. In this article, in addition 
to analyzing the effects of prefiltering and sampling on perfor- 
mance, we seek to find the degradation due to this offset. 

In the advanced receiver being developed for the NASA 
Deep Space Network, a wide range of symbol rates must be 
processed with the maximum available bandwidth and sam- 
pling rate. When the number of samples per symbol is large, 
the loss due to the offset in sampling is negligible. This loss is 
not negligible at high symbol rates, when the number of 
samples per symbol is not large. 

> 

In the digital IDF implementation, we assume that the 
sampling clock is free running at a constant rate, i.e., that the 
sampling clock is not phase locked to a multiple of the symbol 
clock. We also assume perfect symbol synchronization in the 
sense that the receiver has perfect information regarding the 
time that each symbol starts and ends. 

I 

For all of the results presented in this article except Section 
V.E, the symbol period is an integer multiple of the sampling 
period, and the performance is determined as a function of the 
relative phase offset between the sample times and the symbol 
times. In Section V.E, we consider the case for which the sym- 
bol period is not an integer multiple of the sampling period. 
We ignore quantization noise and oscillator instability in this 
article. 

B. Previous Work 

The digital IDF was studied previously by Natali [1] in 
1969. In [1] , a first order low-pass filter was considered and 
the signal-to-noise ratio (SNR) loss was analyzed for different 
time bandwidth products. Only a single symbol was consid- 
ered. Our results match those of Natali for the case of no off- 
set, a single symbol, and using a first order low-pass filter. 

In 1969, Hartman [2] studied the degradation of SNR due 
to intermediate frequency (IF) filtering. He computed the 
SNR degradation due to bandwidth limiting of a biphase mod- 
ulated signal when using an analog IDF. A lower boundary for 
the SNR was derived for different time-bandwidth products. 


In 1976, Turin [3] analyzed the noncoherent digital 
matched filter matched to amplitude-modulated (AM) signals 
in the presence of additive white Gaussian noise and jamming. 
He showed that in the presence of jamming, improvement is 
possible over the analog matched filter by threshold biasing 
and dithering techniques for demodulation. 

In 1978, Lim [4] analyzed the noncoherent digital matched 
filters with multibit quantization, matched to phase-shift- 
keyed signals. 

In 1979, Chie [5] analyzed the digital IDF filter. The per- 
formance of the digital IDF was evaluated in relation to the 
number of bits used by the analog-to-digital (A/D) converter, 
the bandwidth of the prefilter, and the gain loading factor of 
the A/D converter. He also considered quantization error and 
the accumulator length. To study the performance of the digi- 
tal IDF, he considered the symbol error probability resulting 
from hard quantization of the output of the digital IDF. 

It is difficult to determine the exact symbol error probabil- 
ity in the presence of ISI, which is inherent in digital IDF sys- 
tems. The ISI is caused by the low-pass filtering of the input 
signal. Recently, Helstrom [6] and Levy [7] outlined general 
algorithms for approximation of the symbol error probability. 
We should note here that numerous articles and techniques 
address the approximation of the symbol error probability in 
the presence of ISI. Citing all the references on this subject 
is beyond the scope of this article. Furthermore, our intended 
application for this article is the advanced receiver for NASA’s 
Deep Space Network. The receiver output detected symbol 
values are quantized to several bits of accuracy (soft decision), 
as opposed to making hard decisions on the symbols. These 
quantized symbols are then used by the decoder for estimating 
the transmitted bit sequence. In our case, the SNR is therefore 
a more relevant parameter than the symbol error probability. 

We make the fundamental assumption that the sampling 
clock is stable and the receiver symbol clock is synchronized. 
The effects of transmitter receiver clock time-base instability 
on coherent communication systems was analyzed by Chie 
[8] in 1982. The types of time-base instability modeled and 
analyzed are bit jitter and bit jitter rate. 

C. Outline of Article 

In Section II, the IDF is formulated and the average signal 
response and noise variance are derived. In Section III, the 
expression for the loss is formulated. In Section IV, the loss is 
expressed for two special cases: the first order low-pass filter 
and the ideal low-pass filter. In Section V, the numerical 
results are presented for known waveforms and for pseudo- 
random data patterns generated by Monte Carlo simulation. 
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Selections of bandwidth, sampling rate, and asynchronous 
sampling are discussed in this section. This work is summarized 
and conclusions are drawn in Section VI, and Section VII sug- 
gests a direction for future work. 


II. System Description 

The received signal plus noise is denoted by r(f) = s(t - r 0 ) 
+ n(t ), where s(7) is the signal, n(t) is the noise, and r 0 is the 
delay from the transmitter to the receiver. The transmitted sig- 
nal s(t) is a sequence of pulses expressed as 


s(t) =Y,“ k P{t-kT) 


( 1 ) 


At this point we impose no restriction on the shape or dura- 
tion of each pulse p(t ). The input alphabet U is a finite alpha- 
bet with a t e U - {±1, ±2, . . . , ±M} . 

The analog IDF is shown in Fig. 1(a). The analog system is 
an ideal matched filter when p(t ) is a rectangular pulse from 
t = 0 to t = T. It detects the Arth symbol by integrating over 
time kT + r 0 to (k + 1) T + r 0 . The digital IDF is depicted in 
Fig. 1(b). In the digital implementation, an anti-aliasing low- 
pass prefilter is used for filtering the input signal. The filter 
output is sampled, with the ith sample occurring at time 
iT s + r x . The digital IDF detects the kth symbol by summing 
all the samples from t = kT+ t 0 to t = (k + 1 ) T + r 0 . 


We assume that there is perfect symbol synchronization at 
the receiver, so the beginning and end times of each symbol 
are known. For the k th symbol, the “sampling offset” is 
defined by the time difference between the start of the symbol 
and the first sample within the symbol time, i.e., ( iT s + Tj) - 
(kT + r 0 ). The first sample within each symbol time may occur 
anywhere between 0 and T s seconds. A typical symbol wave- 
form and the sampling points are shown in Fig. 2. We seek to 
determine the signal-to-noise ratio (SNR) for the output 
sample of this system. 

Initially, we consider the average signal response of the digi- 
tal IDF; later we consider the noise response. 

A. Average Signal Response 

The response of the low-pass filter to the observed signal 
r(t) is 
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Using Eq. (1) for s(t) we have 
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This signal y(t) is sampled each T s sec at time t = iT s + . 

We denote y(iT s + t x ) as y t . Taking the expectation of Eq. (3) 
conditioned on a given data sequence a , and noting that the 
noise n(t) is assumed to have zero mean, the conditional 
expectation of y { is 

OO ft oo 

E[y. I a] = X) a u I h ( iT s + - f) AS ~ kT " T 0 W 

k- — °° *^-oo 

(4) 

With a change of variable, Eq. (4) can be written as 


E[y t \g] =Yj a k f 

k 


where 6 = t x - r 0 . Let 


h(iT s - kT + 6 - x) p(x)dx 


(5) 


qfk, 6) = f h(iT s - kT+ 6 - x) p(x)dx 

— oo 


( 6 ) 


represent the signal response of the filter at time iT s + t x due 
to a single pulse at time kT + r 0 . For simplicity we denote 
q f (k, 5), as q^k). The total average signal response from 
Eq. (5) for a given fixed 6 may be expressed as 


tfLV/l&fi] =2 a k q f- k ) 


(7) 


Let I k be the set of all i such that the z th sample falls in 
the £th symbol time, i.e., 
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I k = {/: kT<iT s + 5 <(k+ 1)7} 


( 8 ) 



The IDF output for the A;th symbol, denoted by A k , is 


<’> 

• T k 

l€l 

The expectation of A k over the noise, conditioned on a 
and 6 , is 


E[A k \g,8] = EE«. , qfk’) ( 10 ) 

iel k k' 

To further simplify this expression, define the event indi- 
cator function which is 1 if and only if iel k , i.e., 

(l for kT<iT s + 8 <(fc+ l)T 

/,.(5;fc) = J ( in 

( 0 otherwise 

Thus from Eq. (10) we have 


E[A k \g,8] = EE l{8\k)q.(k')a' k ( 12 ) 

i k' 


B. Noise Response 

Next we consider the noise response of the IDF to compute 
the total SNR at the output of the IDF. Let z t denote the 
sampled noise response of the filter at time iT s + r t . 

2 ,. = f n(£)h(iT' + r x -m ( 13 ) 

oo 

Since the IDF is a linear system, the variance of .4*. con- 
ditioned on a is equal to the variance of the response of the 
kt h symbol due to noise alone, i.e., it is independent of s(t). 
The variance of A k is 


E[zz.] = R z (i-0 = y J 


( 15 ) 


where R z is the autocorrelation of 


III. Definition of Signal-to-Noise Ratio Loss 

The analog IDF of Fig. la is the optimum matched when 
p(t) = 1 for 0 < t < T and zero otherwise. The SNR is defined 
at the IDF output as the ratio of the square of the mean to the 
variance. Denoting SNR^ for the analog IDF, it is well known 
[ 10 ] that 


SNR ^ 


24 2 r 

N n 


(16) 


We assume with no loss of generality that the signal ampli- 
tude A = 1 , so SNR^ = 2T/N 0 . 

The SNR at the output of the digital IDF is denoted by 
SNR d . We compare the performance of the digital IDF with 
the analog IDF by considering the ratio 


a SNRg 

7 “ SNR* 


(17) 


We define SNR D at the output of the IDF as the ratio of 
the square of the conditional mean to the conditional variance 
of the output, conditioned on a given sequence and offset 
value. From this definition of SNR, 


{E[A k \q,8]f 
Sm D ~ var^la, 5] 


Then we have 


7 = 


N o 

2T S NRd 


(18) 


vai[A k \ a, 6 ] = EE 1.(8; k) 1.(8, k)E[zz.] 


Note that this variance does depend on 6 and k , because the 
number of samples occurring in the k th symbol depends on 5. 
Using Eq. (3) and noting that E[n(i)n(r)] = NJ 2 5 0 (f - r), 
where here 6 0 (*) is the Dirac delta function, we have 


In the remaining sections, 7 dB = 10 log 10 ( 7 ) (dB) is com- 
puted for various filters and data patterns. Normally 7^1, 
because the digital IDF has a loss with respect to the analog 
IDF. The degradation or loss in decibels is the negative of 7 dB , 
and minimum loss corresponds to the maximum attainable 7 , 
which typically approaches one ( 7 dB = 0 dB). Maximum loss 
is unbounded. In some cases, for a given data pattern a, there 
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is a gain in SNR (e.g., all-ones sequence), in which cases y > 1 
and the degradation is negative in decibels. 


IV. SNR Performance for Special Cases 

In general, the pulse shape p(t) may be chosen to take 
numerous shapes (e.g., raised root-cosine). In some cases, 
it is chosen to extend over more than one symbol duration, 
e.g., for partial response signaling (sometimes referred as 
correlated coding or controlled intersymbol interference). 
For bandwidth-limited channels, the pulse shape and duration 
are selected to increase the bandwidth efficiency of the com- 
munication system. 


and W is the radian cutoff frequency of the filter. The filter in 
this case is causal and therefore physically realizable, i.e., 
h(t) = 0 for t < 0. Using Eq. (20) to evaluate qfk) from 
Eq. (12), 


(1 _ e -fcW s - fcr) ) . kT<iT+5 <(k+ l)T 


qfk) = < e 




; (*+ \)T<iT +6 
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We consider only non-overlapping rectangular pulses 
throughout the rest of this article, since this pulse shape has 
traditionally been used for NASA’s deep space missions. It is 
pointed out that the results of the previous sections, namely 
the expressions for average signal response in Eq. (12) and the 
noise variance in Eq. (15), hold regardless of the pulse shape or 
its duration. 

In the case of the rectangular pulse we simply have 


pit) = 


then from Eq. (6) 


1 0 <t<T 


0 otherwise 


Qfk) = / KiT s -kT+ h-x)dx (19) 


To derive the equation for the noise variance, we use the 
expression in Eq. (15), resulting in 
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To obtain the relative performance of the digital IDF when 
using a first order low-pass filter, we use the signal response ex- 
pression in Eq. (21) and the noise variance in Eq. (22). Then, 
from Eq. (18), y is 
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E[A k \a,f>] = ££ lfb-k)qfk')a k 


where qfkT) is defined in Eq. (21). For a known data pattern, 
Eq. (21) may be evaluated for different k , 5, and time band- 
width products {WT). 


It is useful to write Eq. (19) as 


r(k+ 1U 
J kT 


h(iT +b-x)dx (20) 


We now consider two different low-pass filters, one causal 
and one non-causal. 


B. Ideal Low Pass Filter (Non-Causal) 

The ideal low-pass filter with unit gain and low-pass band- 
width W Hz is non-causal with the impulse response 


h(t) = 2 W = 2W sin c(2nWt) (24) 

Z7T W t 


The expression for the signal response in Eq. (20) becomes 


A. First Order Low-Pass Filter (Causal) 

The impulse response for a first order low-pass filter is 
h(t) = We~ Wt u(t ), where u(t) denotes the unit step response 




j r(k+i)T sin 2v + 5 - *) 


0T+8-*) 


dx (25) 
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It is possible to express Eq. (25) in terms of 
, C x sinu , 

= J' — * 

as 

qfk) [Si(2irW(iT-(k+l)T+6)) 


filter such as that of Butterworth or Chebyshev [11] and [12] 
does not greatly influence the results, since the realizable filter 
can be considered an approximation to the unrealizable filter 
with finite group delay [12] . 

All the computational results were obtained by explicitly 
evaluating Eq. (30) and Eq. (23) for arbitrary input signal 
sequences a. The Si(') function in Eq. (30) was implemented 
using a 500 point look-up table provided in [9] . 


-Si(2irWQT -kT + S))] (26) 

To find the noise variance, it suffices to note that the noise 
spectral density at the output of the filter is 

for all | / | < W 

(27) 

otherwise 

and thus the autocorrelation function is 



R n (r) = N q W 


sin 27 t W t 
2n W r 


(28) 


The simulation result for the first order low-pass filter SNR 
[Eq. (23)] is not included in this article. It was, however, com- 
puted by the authors to verify the results against Natali’s re- 
sults [1] , as was pointed out in the introduction. 

For Figs. 4 through 13, the letters a through / correspond 
to the following time bandwidth products (WT): a = 2, 
b = 1.75, c = 1.50, d = 1.0, e = 0.75, and /= 0.5, unless stated 
otherwise. 

A. Spectral Analysis 

To gain a better understanding of the digital IDF, we inves- 
tigated the spectral properties of the signals processed in the 
IDF. We considered two cases. For the first case, the data sig- 
nal s(t) is the alternating +1, -1 sequence with the Fourier 
series 


Thus, the noise variance at the output of the IDF can be 
expressed from Eq. (28) as 




var[.4 fc l6] = tEE sinc(2H'(i'-./)2^) (29) 


. T k . r k 
tel jel 


sit) = £ £ cos 2nf 0 (2i+ 1 )t (31) 

1=0 

where f Q = 1/2T. We refer to this case as case 1. Note that the 
signal has energy only at odd harmonics of^. 


The performance y can be evaluated from Eq. (18), which 
yields 



For case 1 and the ideal filter, the bandwidth W determines 
the number of odd harmonics that pass through the filter. 
When the time bandwidth product WT is greater than 1 .5 and 
less than 2.5, the first two odd harmonics pass through the fil- 
ter; and when 0.5 < WT <1.5, only the first harmonic is 
passed through the filter. 

The second case, case 2, is the binary random waveform 
with rectangular pulse shape, which has the autocorrelation 
function [13] R s (r): 


V. Numerical Results 

| 


|r| <T 


Numerical results are presented only for the ideal low-pass 
filter. The motivation to consider the ideal low-pass filter is 
to eliminate aliasing in an ideal manner. The use of a realizable 

= 

to 

otherwise 

(32) 
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with power spectral density 


S s (f) = T 


sin 2 (vfT) 
(rr/T) 2 


(33) 


The noise process n(t) and the transmitted signal s(t) are 
mutually independent. The autocorrelation function of the 
observed signal r(t) is 


R r (r) =R s (r)+R n (T) (34) 


The spectral shapes of the transmitted signal s(f) and the 
noise process n(t) are shown for case 2 in Fig. 3. The time- 
bandwidth product WT determines the number of lobes of 
the signal spectrum Eq. (33) passed through the filter. 

B. Digital IDF Output Noise Variance 

The noise variance at the output of the IDF for the case 
when there are N samples in the symbol is 

N- 1 

varM fc | a,S] = NR n (0) + £ 2(N - i) R n (iT) 


This variance reduces the NR n ( 0) when the samples are 
independent. This is possible by choosing the sampling period 
T s = 1/2 W, which is evident from Eq. (28). 


C. Performance Loss Versus Offset and 
Bandwidth for Known Waveforms 

In this section, we consider known signals as input to the 
digital IDF. The output of the ideal low-pass filter depends on 
both past and future inputs. To approximate this, we consider 
21 -symbol blocks, and the 11th symbol is analyzed for each 
data pattern. A block of 21 symbols was found to be suffi- 
ciently long to analyze the IDF for different data patterns. 
This is reaffirmed in the following section (V.C.l) by consider- 
ing the spectrum of the sampled waveform for the alternating 
data pattern. Five different data patterns are considered: 

(1) Alternating data pattern: a = (-1, +1, -1, +1 , -1, . . .) 

(2) Single pulse: a = (-1, -1, -1, . . ., -1, +1 , -1, . . . ) 

(3) Two ones: a = (-1, +1, -1 , +1 , . . ■ , +1, -1 , +1 , +1 , -1, 
+1,...) 

(4) Three ones: a = (+1,-1,+1,-1,...,+1,-1,+1,+1, 
+ 1,-1, +1,-1,. . .) 

(5) All ones: £ = (+1, +1, +1, ... ). 


We analyze the waveform and SNR for each pattern for dif- 
ferent WT, particularly for when there are four samples per 
symbol, i.e., T s = T\ 4. 


In the following sections and related figures, the offset 6 is 
defined as the length of time from the start of the 11th sym- 
bol to the time when the first sample of the 11th symbol 
occurs. 

1. Pattern 1, alternating data. The sampled waveform for 
the 11th symbol, a -1 pulse, is shown in Fig. 4(a) for WT - 2. 
In Fig. 4(a), for every sampling offset value, with increments 
0.05 T s , a unique English letter (a through t) is used to indicate 
the point at which the sample occurs. Every letter occurs four 
times, corresponding to the four samples per symbol. 

In Fig. 4(b), the anti-aliasing filter output waveform is 
shown for different time bandwidth products. The sampled 
waveform approaches a sinusoid when 0.5 < WT <1.5, and 
approaches the sum of a sinusoid and its third harmonic when 
1.5 < WT < 2.5. This indicates that the sequence of 21 data 
symbols is sufficiently long to approximate the infinite alter- 
nating sequence, since the waveforms agree with harmonic 
properties discussed earlier in Section V.l. 

The degradation is depicted in Fig. 4(c) for different offset 
values, when WT ranges from 0.75 to 2.0. For all WT cases, the 
worst case occurs when 5=0, with the value of the worst case 
ranging from 1.6 dB at WT = 0.75 to 1.96 dB at WT - 1.75. 
To find out how much loss is due to sampling, the loss for 
T/T s = 20 is shown in Fig. 4(d). As expected', since the sam- 
pling rate is high, the loss depends mainly on WT , not on 5. 
For the case when WT - 1, the worst case loss is decreased 
from 1.92 dB at T s = 27 4 to 1.10 dB for 7; = 77 20. When T s = 
7720, the digital IDF almost approximates the analog IDF 
(with finite bandwidth). Hence, it can be deduced that about 
0.82 dB loss results from sampling with T s = 77 4, and 1.1 dB 
loss results from the bandwidth limiting of the received signal. 

2. Pattern 2, single pulse. We consider here the second 
pattern and compare it to the first pattern. The results are 
compared to those of Hartman [2] , who considered the same 
pattern. 

Hartman points out that when the time-bandwidth product 
is an even integer, the maximum loss results when the binary 
waveform is a single 1 preceded and followed by all -Ts (or 
vice versa). He also shows that when the time-bandwidth 
product is an odd integer the maximum loss results when 
the binary waveform consists of alternating Ts and -l’s, 
respectively. 
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In Fig. 5(a) the resulting waveform for this sequence is 
shown for different WT values. In Fig. 5(b) the degradation is 


shown as a function of the offset 5. Performances of the alter- 
nating data pattern and a single pulse are compared in Fig. 6, 
for WT = 2. We find that the loss is virtually the same for 
the two sequences, and the alternating sequence results in 
slightly more degradation, approximately 0.01 dB for the same 
values of offset. This is in slight disagreement with Hartman 
[2] . For WT = 1, the worst case loss is 1.77 dB for the single 
pulse, and 1.92 dB for the alternating pattern, which agrees 
with Hartman [2] . 

3. Pattern 3, an asymmetric pattern. The third data pat- 
tern is (-1, +1, . . +1, -1, +1, +1, -1, +1, . . .) considered. 

The filtered output waveforms shown in Fig. 7(a) and Fig. 
7(b) depict the corresponding degradation. Since the wave- 
form is asymmetrical about the Tj 2 point, the loss is also 
asymmetrical about the TJ2 point. 

4. Pattern 4, alternating except for three ones. The loss 
is shown for this pattern in Fig. 8. The sampled waveform is 
not shown in this case, since it is almost constant throughout 
the symbol time. The performance curve in this case is sym- 
metrical, since the symbol waveform is symmetrical around 
the 772 point. In this case, the worst case loss does not occur 
for 5 = 0 or 5 = T s , which is different from the previous cases. 
Also, for cases b, c, and d (WT - 1.75, 1.5, and 1.0), there is a 
gain rather than a loss in SNR for all or most offsets. This 
occurs because the ISI happens to aid in these cases. 

5. Pattern 5, all ones. Cases 3 and 4 result in a slight gain 
in the SNR for certain offsets and time bandwidth products. 
To investigate this point, the loss curve in Fig. 9 is shown, when 
a step function (all 1 ’s pattern) is used as the symbol sequence. 
In this case when WT < 1 is selected, the performance curve 
indicates a constant gain. This is a consequence of the expres- 
sion for 7 in Eq. (30). By decreasing the anti-aliasing filter 
bandwidth to WT < 1, the amount of filtered noise power 
decreases. The signal power does not significantly decrease for 
long periods of constant data. This results in a performance 
gain for these cases. 

D. Selection of Bandwidth and Sampling Rate 

In implementation of the digital IDF, for a given sym- 
bol rate T y a reasonable criterion is to select the sampling 
period T s and the bandwidth W such that the average loss 
is minimized, where the average is over the ensemble of all 
possible data patterns. One could consider either averaging the 
loss over all offsets or the worst case offset. Both bandwidth 
and sampling may also be restricted by hardware or other 
considerations. 

It is well known, when applying the orthonormal set of 
radial prolate spheroidal functions, to expand a band-limited 


signal in the function space of finite energy signals, that the 
number of orthogonal dimensions (eigenvalues) necessary to 
describe a band-limited function over a T-sec time interval is 
approximately 2WT [12] . 

The Nyquist sampling theorem [11] requires that the sam- 
pling period satisfy the inequality 


or, equivalently, the number of samples/symbol must satisfy 
the inequality 

J> 2WT (37) 

S 

The digital IDF asymptotically approaches the analog IDF 
when W y and T/T v approach infinity. This leads to an infinite 
time-bandwidth product which is unrealizable in a physical 
system. 

The bandwidth W should be selected such that most of the 
signal power is passed through the filter. Furthermore, it is 
clear from Eq. (36) that as the bandwidth is increased, the 
required sampling period decreases, and more samples per sym- 
bol are required. 

To select the best filter bandwidth given a symbol rate, one 
must also consider the current hardware constraints of tech- 
nology. The two most important constraints are the speed and 
accuracy of available A/D converters, and the speed and accu- 
racy of signal processing hardware such as multipliers and 
accumulators. Therefore, the sampling rate and the filter band- 
width must be selected such that the hardware is practically 
implementable. 

1. Optimum choice for T s and W. Given a symbol rate, 
the necessary steps for the design engineer are: 

(1) To determine the maximum practical limit to the 
sampling rate or, equivalently, the maximum number 
of samples that is possible during each symbol time. 

(2) Once the sampling period T s is determined, then the 
bandwidth W must be selected such that the inequal- 
ity [Eq. (37)] is satisfied, preferably with equality. 
This results in passing maximum signal power into the 
digital IDF, while satisfying the sampling criterion. 

Investigating the SNR loss when Eq. (37) is satisfied with 
strict equality is the subject of the study in the following 
section. 
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2. Monte Carlo simulation. We consider the class of sym- 
bol patterns that are random binary sequences selected from 
an equally probable binary alphabet. In Figs. 10 through 
13, the average loss for binary vectors of length 4620 bits, 
generated from a binary symmetric source, is computed for 
various time-bandwidth products and sampling rates. 

The averages are computed by breaking the 4620 bit vector 
into 220 blocks of 21 symbols each and computing the loss of 
the 11th symbol for each block. It was confirmed previously 
that 22 symbols are sufficiently long to approximate an infi- 
nite length sequence, for simulation purposes, when the 
11th symbol is being analyzed. 

The stopping rule for computing the average loss is to repeat 
computing the loss for every block until the accumulated loss 
does not change by more than 10~ 4 for 10 consecutive blocks. 
At this point the average loss is computed by dividing the 
accumulated loss by the number of blocks that have been 
processed. This approach has the advantage of averaging over 
the ensemble of ISI patterns. Furthermore, the incremental 
loss computed for a given block of data patterns is indepen- 
dent of the loss computed for the previous data pattern. 

The results confirm that the minimum loss averaged over 
the ensemble of all possible test patterns is achieved when 5, 
the sampling offset, is TJ 2. This conclusion is only true if the 
sampling rate is an integer multiple of the symbol rate. This is 
clarified when we consider the case in which the symbol per- 
iod is not an integer multiple of the sampling period. 

Figure 10(a) depicts the average loss versus offset for T s = 
7 1 / 2 and for several different WT products. For WT = 1, the 
worst case SNR loss in this case is about 2.4 dB, at 5 = 0. 
When WT > 1 the inequality of Eq. (37) is violated, because 
the sampling rate is too low for the bandwidth, and aliasing 
occurs. The curves for WT > 1 are depicted to exhibit the loss 
due to aliasing. Only curves d and e , with WT = 1 and 0.75, 
satisfy Eq. (37). For the best offset, 6 = 0.5, which corre- 
sponds to phase locked sampling, the loss is only 0.7 dB. 

Figure 10(b) depicts the case when T/T s = 4. In this case 
Eq. (37) is satisfied for all WT considered. For most values of 
the offset, the loss decreases uniformly as the time bandwidth 
product approaches 2. Thus, selecting W = 2/7 is the most 
appropriate choice, which satisfies Eq. (37) with strict equal- 
ity. The worst case loss occurs when 6 = 0, and the corre- 
sponding loss is approximately 1.2 dB for all WT cases. The 
widest bandwidth is recommended because it minimizes the 
average loss, averaged over all offsets. 

Figure 10(c) depicts the case when T/T s = 8. In this case the 
minimum loss (at 8 - T s /2) is virtually unchanged from the 


case when T/T s = 4. This loss is approximately 0.3 dB to 0.4 dB 
for WT - 1.5 and 2.0. For the widest bandwidths, WT - 1.5 
and 2.0, the worst case loss is improved by 0.6 dB, from 1 .2 dB 
to 0.6 dB, over the case when T s = 7/ 4. This improvement 
results from doubling the number of samples per symbol 
at BT=2. 

The conclusion is that sampling at twice the Nyquist rate 
reduces the worst case loss by 0.6 dB compared to sampling 
at the Nyquist rate, for WT - 2. Sampling 16 times the symbol 
rate with WT - 4 [see Fig. 1 1(a)] results in the worst case loss 
of 0.3 dB, of which 0.2 dB is due to bandwidth limiting and 
only 0.1 dB is due to sampling and the IDF. 

Two cases are considered to investigate the case when WT is 
selected to be greater than two. In Fig. 1 1(a) with WT - 4, and 
in Fig. 11(b) with WT = 8, the average degradation is shown 
for different T/T s ratios. As the number of samples/symbol is 
increased in Fig. 11(a), the variation in the loss due to the 
offset becomes negligible and both the worst case and best 
case losses approach approximately 0.2 dB. This indicates that 
the degradation due to the offset becomes negligible as the 
number of samples per symbol is increased. 

Regressing to the practical scenario, suppose that given 
practical hardware constraints, it is possible to obtain a maxi- 
mum of eight samples per symbol. To find the best choice for 
the bandwidth, we have compared the loss for WT = 2 and 4 in 
Fig. 12. It is apparent that in the worst case offset (6 = 0), 
both cases result in equal loss of 0.6 dB, but the minimum 
degradation is better by 0.2 dB when WT = 4. Thus, this 
choice of bandwidth is better when T s = 77 8. 

3. Summary of Monte Carlo simulation results. To sum- 
marize the results for simulation, it is concluded that 0.75 < 
WT < 1 is appropriate when T s - T/2. When T s = 774 it is 
recommended to use the widest filter that eliminates aliasing, 
i.e., W = 2/7". Then sampling is at exactly the Nyquist rate. 
This is better than selecting WT < 2 and sampling higher than 
the Nyquist rate. When eight samples per symbol are available, 
i.e., T s = 77 8, it was deduced that the worst case loss is the 
same for W = 2/7 and W = 4/7. However, the minimum aver- 
age loss is lower for wider bandwidth, W = 4/7. When the 
sampling rate is limited, the widest bandwidth that satisfies 
Eq. (37) should be used. Sampling at twice the Nyquist rate 
is recommended when the bandwidth is limited. 

E. Effects of Asynchronous Sampling 

Normally, in practical systems the sampling clock is not an 
exact multiple of the symbol clock. This occurs when the sam- 
pling is at a fixed rate and the symbol rate varies, for example, 
with Doppler shift. 
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Figure 1 3 depicts the average loss for a random data pattern 
when T/T s = 4.5. Some symbols have four samples, and some 
have five. When the number of samples is four (or when the 
offset 6 < TJ 2), the degradation is approximately G.3 dB less 
than when the number of samples is five. This is true for all 
WT considered. The average loss for five samples per symbol is 
due to additional noise power in the last sample, since the last 
sample occurs near the transition point of the filtered symbol 
waveform. Despite this variation in loss with the number of 
samples, the worst case loss for T/T s = 4.5 is only 0.7 dB for 
WT - 1.5 and 2.0, which is less than the worst case loss of 
1 .2 dB when T/T s = 4. Even though there is a variation in loss 
with number of samples per symbol, a higher sampling helps 
performance rather than hurting it. 


VI. Summary and Conclusions 

The performance of the digital IDF was studied in this 
article, with special attention to the effect of offset in sam- 
pling. We derived the expressions for the signal response, noise 
variance, and IDF output SNR for general pulse shapes and 
anti-aliasing filters. These results were then specialized to the 
rectangular pulse shape and the ideal low-pass filter. Perfor- 
mance was evaluated in terms of loss in SNR compared to an 
analog IDF with infinite bandwidth. 

Performance was evaluated for five different known data 
sequences, and for random sequences. Observations regarding 
the worst case degradation and the effects of the offset in the 
sample time were made. It was concluded that when the sys- 
tem is bandlimited, i.e., the time bandwidth product WT is 
finite, a sampling rate in excess of the Nyquist rate (2 BO 
should be chosen. For example, when the bandwidth is only 
twice the symbol rate, sampling at twice the Nyquist rate 
results in a 0.6 dB improvement over sampling at the Nyquist 


rate, when the worst case loss is the performance criterion. 
When the sampling rate is limited, i.e., T/T s is fixed the band- 
width W should be selected such that W = 1/2 

The above results are for the case when the sampling rate is 
not phase locked to the symbol rate. For an integer number of 
samples per symbol, it was shown that the degradation due to 
offset in sampling is minimized when the offset is half the sam- 
pling period, i.e., 5 = TJ2. For four samples per symbol, the 
worst case loss is 0.9 dB greater than the best case. This means 
that phase locked sampling is 0.9 dB better than the worst case 
for the non-phase-locked sampling, when the ratio of symbol 
rate to the sampling rate is small (<4). 

The loss due to the offset becomes negligible when the 
number of samples/symbol becomes sufficiently large. With 
WT > 4, T s < 1/2 W, the digital IDF always performs within 
0.6 dB of the analog IDF, for the random data pattern and the 
worst case sampling offset. 

The effect of a non-integer ratio of sampling rate to symbol 
rate was also studied, for the case when T/T s = 4.5. For a given 
time bandwidth product, the worst case loss is lower than the 
case when T/T s = 4. Thus letting the sampling rate be a non- 
integer rate relative to symbol rate did not degrade the 
performance. 

VII. Direction for Future Research 

For the Deep Space Network receiver, the losses incurred 
for the digital IDF at high symbol rates are undesirable. To 
reduce the loss, it is possible to weight each sample in the IDF 
such that the loss is minimized. That is, instead of performing 
a simple summation, a weighted summation is performed. This 
will be considered in a later article. 
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Glossary of Terms 


T 

s 

- Sampling time in seconds 

m 

- The transmitted signal 

T 

- Symbol time in seconds 

t) 

- Additive white Gaussian noise with flat spectral 

W 

- Filter bandwidth in hertz 

•it) 

density NJ2 
- The received signal 

a. 

i 

- Transmitted symbol 

y{t ) 

- The output of the low-pass prefilter 

P(0 

- Pulse shaping waveform 

m 

- The transfer function of the low-pass prefilter 


- The sampled output of the Integrate-and-Dump Fil- 

T o 

- Transport lag from transmitter to receiver in seconds 


ter at time k 

T 1 

- Sampling delay in seconds 

y* 

- The sampled output of the pre-filter 


- Autocorrelation function of signal x(t) 
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Fig. 1. Analog and digital Integrate and Dump Filter (IDF): 
(a) analog IDF; and (b) digital IDF 



Fig. 2. Offset in sampling 



Fig. 3. Signal and noise spectral density for binary random 
waveform: (a) spectral density of transmitted signal; (b) spectral 
density of the noise; and (c) frequency response, the ideal filter 
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Fig. 6. Degradation comparison vs. offset for alternating and 
pulse data pattern 
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